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Introduction 

Aerobraking has been proposed as an efficient means of decelerating spacecraft for planetary 

missions. 1 Most current aerobrake designs feature a blunt forebody shielding the payload from 
the intense heat generated during atmospheric entry. Although this forebody will absorb the larg- 
est portion of the heat pulse, accurate prediction of heating in the near wake is of great impor- 

tance, since large local heating values can occur at points of shear-layer impingement. 

In order to address the various issues associated with these blunt-body wake flowfields, the 
Advisory Group for Aerospace Research and Development (AGARD) formed Working Group 1 8 
in 1992. One of the objectives of this activity was to examine real-gas effects in high-speed flow- 
fields around a 70° blunted cone; the primary dimensions of this geometry are shown in Fig. 1. To 
date, many researchers have conducted experiments using this geometry in various facilities, such 

as the Large Energy National Shock (LENS) tunnel at Cubric/Calspan 3 and the HEG shock tunnel 
at DLR-Gottingen. 4 Several computational studies have also been conducted in concert with these 
tests. 5-6 

Many of the experimental results have indicated the possible presence of a transitional shear 
layer through a large increase in heat transfer downstream of the reattachment point. The presence 
of transition could in fact lead to much higher peak heating than if the separated flow is entirely 

laminar or turbulent. 7 In the shock-tunnel tests, however, it is difficult to separate such viscous- 
flow phenomena from real-gas effects. In order to help make this distinction, Horvath et al. 
recently conducted a set of experiments in the NASA Langley 20-Inch Mach 6 Tunnel, and com- 
pared the results to laminar Navier-Stokes calculations. 8 They found heat-transfer distributions 
similar to those obtained in the high-enthalpy facilities, with the measured peak heating along the 
sting support markedly greater than that predicted by the laminar computations. These trends 
point to the need to find transitional and turbulent computational solutions for these flowfields. 
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The objective of this work is to assess the requirements for transition and turbulence modeling 
in blunt-body wake flows. 


Present Approach 

Our approach relies on the implementation of appropriate transition and turbulence closure 

models in the nonequilibrium flow solver developed by Olynick. 0 This is a 5-species, 3-tempera- 
ture algorithm for the full Navier-Stokes equations. It utilizes Roe’s invscid flux-difference 

splitting 10 with variable extrapolation for high-order extension, and the LU-SGS diagonal 


implicit algorithm of Yoon 1 1 for time integration. An attractive feature of this implicit method is 
that the computational expense associated with the method increases only linearly with the num- 
ber of partial differential equations being solved. 

In order to include the influence of turbulent fluctuations in the algorithm, the conservation 
equations are averaged using the standard Favre-averaging technique, giving rise to additional 
unkowns commonly referred to as the Reynolds stresses and Reynolds heat flux. These terms are 
modeled using Boussinesq’s approximation, so that 
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where the overbars indicate time-averaged quantities, the tildes denote Favre-averaged quantities, 
the double primes indicate fluctuations, and k = is the turbulent kinetic energy. The tur- 


bulent transport properties are assumed to be functions of turbulent length scales, as well as the 
turbulent kinetic energy. Separate length scales are used to define the turbulent viscosity and ther- 
mal conductivity, leading to a variable turbulent Prandtl number. Evolution of the turbulent 
kinetic energy is presently governed by a single additional partial differential equation, discussed 
in Ref. 12. 

This one-equation model is augmented by a high-speed transition model due to Warren et 

a/. 13 The transition model incorporates first-mode disturbance scales extended to compressible 
flows using the reference temperature method, as well as second-mode disturbance scales often 
present in supersonic flowfields. 


Results 

Preliminary results for a freestream Mach number of 6 and a freestream Reynolds number of 2 
million have been obtained on the 125 x 90 grid shown in Fig. 2; because of the relatively low 
freestream Mach number and the desire to separate viscous effects from real-gas effects, chemical 
reactions are disabled for these computations. Figs. 3 and 4 compare numerical heat-transfer 
results for both laminar and transitional solutions to experimental data from Ref. 8 for this flow 
condition. Note that the transition point was chosen based on the location of the increase in the 
experimental heat-transfer rate. The forebody results in Fig. 3 show very good agreement with the 
experimental data. However, in Fig. 4, it is evident that there is a substantial discrepancy between 
the predicted heating values and the experimental measurements. Most importantly, the transi- 
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tional calculation does not show nearly as large a heating rise as is seen in the experimental data. 
All available transition models, including that developed in Ref. 13, address transition in one 
shear layer and all require a specification of a transition onset location. The wake flow of a body 
mounted on a sting is characterized by the interaction of a free shear layer and the sting boundary 
layer. Based on the experimental data presented in Ref. 8, transition is taking place within the 
shear layer. Thus, specifying a transition point on the sting is not going to produce an accurate 
description of the flowfield in the wake region. Therefore, existing transition models are incapable 
of handling such flows. 

Future Work 

We plan the implementation of a newly developed transition model 14 that determines transi- 
tion onset in the flow as part of the solution. Further, the model is complete in the sense that no 
length scales need to be specified. All it requires are initial and boundary conditions. 
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WFigures 


70° Blunted Cone with Sting 

R b =7.62 cm, RyR^O.50, R,/R b =0.25, \-JR b =6 

Location s/R n X, cm Y, cm 

1 0.349 0.230 1.303 

2 2.043 2.438 7.369 

3 2.166 2.796 7.620 

4 2.323 3.177 7.239 

5 3.557 3.177 2.537 

6 3.817 3.809 1.905 

7 15.652 48.897 1.905 
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Figure 1. Blunted-cone geometry 




Figure 2. Grid for Navier-Stokes calculations 
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